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Abstract 

We study the resummation of large QCD radiative corrections up to the next to lead- 
ing logarithmic accuracy to the process 77 — s- bh\ i.e., we resum logarithms of the type 
a^glr?^ m? / s and af ln^^~^ m^/s (m is the quark mass). The only source of all the log- 
arithms to this accuracy is the off-shell Sudakov form factor included into the triangle 
topologies of the one-loop box diagram. We prove that any other configurations of dia- 
grams to this accuracy, either cancel in subgroups or develop a universal on-shell Sudakov 
exponent due to the final quark anti-quark lines. We study the mechanism of cancella- 
tions between the different diagrams, which leads to the simple resummed results. We 
show the cancellation explicitly at three loops for the leading and at two loops for the 
next-to-leading logarithms. We also point out the general mechanism responsible for it, 
and discuss how it can be extended to higher orders. 
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1 Introduction 



Future Linear Colliders are expected to reveal the answers to many questions of modern particle 
physics. One of these concerns the the physics of the Higgs particle and origin of electroweak 
symmetry breaking. The neutral scalar Higgs boson is an important ingredient of the Standard 
Model (SM) and is the only SM elementary particle which has not been detected so far (see 
for a review 0). The lower limit on rriH of approximately 113.5 GeV at 95% c.l. has 
been obtained from direct searches at LEP [Q. Current experiments are concentrating on the 
possibility of finding a Higgs particle in the intermediate mass region 113.5 < ttih < 150 GeV. 
In this region it decays mainly to a bb pair. 

The photon mode of the future Linear Collider (LC), namely the collisions of the energetic 
polarized Compton photons, will be used for the production and for the study of the Higgs 
particle. In the intermediate mass range, the main production process is 

77 i7 ^ bb. 

QCD as well as electroweak radiative corrections to this process have been studied very well 
and have been found to be small in this region The main challenge, however, is to get under 
control the background process 

77 — > bb, 

which gets extremely large QCD corrections. 

In this paper we discuss the process of the quark anti-quark production in the photon 
mode of the LC, 77 bb. The amplitude for this process in the scalar channel contains large 
double logarithms at |m| ^ m^. At very high energies, the large logarithms spoil the 

perturbative predictions. Therefore it is mandatory to develop a clear resummation procedure 
for this double and single logarithmic terms. Let us stress, that the main interest in the 
process 77 bb comes from the fact (but is not limited to it) that it represents the dominant 
background for the production of the Higgs particle, •yy ^ H ^ bb [|, 0|. In fact, our 
motivations for this study are twofold: 

(1) A detailed study of the process 77 bb is very important due to phenomenological reasons 
mentioned above; and 

(2) in addition to that, the quark anti-quark production in photon collisions , being one of the 
simplest process in QCD, is important in its own right, e.g. for studying and understanding 
QCD effects. 

The Born cross section for the polarized bb production in the scalar channel (at 7^ = 0), 
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where the Higgs will be studied as well, is suppressed by ^ H) > (here s is the center of mass 
energy of the initial photons). However, the perturbative QCD corrections contains the large 

double logarithms of the form p = asln^(^), which give a contribution to the cross section 
which is of the same order as the Born contribution at high energies. The presence of the large 
correction was noticed by Jikia in 0. The double logarithmic nature and the origin of these 
corrections were studied in M. The authors studied the process to one and two loop accuracy. 



Later, the form of the resummed results for the double logarithms has been argued in |[10 
These authors also claimed that the double logarithms have a"non-Sudakov" origin. 

In this paper, we first present an alternative way of understanding the resummation proce- 
dure for the double logarithms. The general idea of our approach, is that the only source of 
double logarithms is the off-shell Sudakov form factor included in the triangle topologies of the 
one-loop box diagram. We have proved that the other types of the higher loop diagrams will 
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either cancel in subgroups or develop a universal on-shell Sudakov exponent due to the final 
quark anti-quark lines. In addition to that, (1) we extend this analysis to the next-to- leading- 
logarithmic (NLL) accuracy, and (2) we study the mechanism of cancellations alluded to above, 
between the different diagrams which leads to very simple resummed results. We demonstrate 
the cancellation explicitly at three loops and point out the mechanism responsible which then 
allows for the generalization to higher loops, all up to the next to leading logarithmic accuracy. 
As an aside, we argue that all the large logarithms up to the next to leading level are related 
to the Sudakov ones. This includes not only the leading ones of the form a^^lr?^ m? / s but 
also the next to leading ones of the form a^ln^^~^m^/s (m is the quark mass). It is this fact 
together with an understanding of the cancellation mechanism which allows us to develop an 
easy resummation procedure. 

The paper is organized as follows. In section 2.1 we discuss the one loop diagrams and 
explain which topologies are important. In section 2.2 we present a procedure for the resum- 
mation of the double logarithms. In section 2.3 we extend our analysis to the next-to-leading 
logarithmic level. In this section we consider all possible topologies and give the final result 
of the resummation. Section 3 is devoted to the study of the different cancellations which are 
responsible for the simple resummed results of section 2.2 for the double logarithmic case. In 
section 4 we tackle the case of the form factors to the next to leading logarithmic accuracy and 
justify the results given in section 2.3. The paper ends with some conclusions and discussions. 

2 The resummation up to next to leading logarithms 
2.1 One loop results 

We study the process of quark anti-quark pair production in photon-photon collisions 

+ l{k2) ^ q{pi) + q{p2) (1) 

in the scalar channel Jz = 0, at very high energies compared to the mass of the quark m, and 
at large scattering angles 

|s| ~ \t\ ~ \u\ > m^. (2) 

The radiative corrections to this helicity amplitude contain large double and single QCD loga- 
rithms of the form ln^(s/m^),p = 2, 1, in the limit eg . ([T7|JT8| ) . 

It has been shown in that the large double logarithms (DL) have a Sudakov-like nature 
and can be extracted from the total result by identifying special kinematic regions. It was 
demonstrated that, only the box diagram contributes, and that this box diagram can be reduced 
to three different effective diagrams with triangle topology. We label them as topologies A,B,C 
following 0, see Fig.(|l|). These effective diagrams result from the box diagram when one of 
the four propagators is hard, two collinear and one soft ( the soft line can be either the gluon 
or the quark propagator). Hence one can shrink the hard propagator of the box diagram to 
a point, because it does not depend on the soft loop momenta anymore. The hard momenta 
would be just an ultra violet cutoff in this situation. 

In order to introduce notations and to get an idea about the structure of the DL, let us 
consider one loop calculations in the DL approximation. Also, we will use these results in the 
next section for the construction of the DL-resummed result. 
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A B Ci C2 

Figure 1: Triangle topologies obtained from the one loop box diagram. 



As is well known, the double logarithms come from the region of soft loop momenta. Hence, 
one can always write a factorized form for the amplitude in the DL approximation: 

Ma = Mgorn • Fa, Mb = Mg^^^ ■ F^, Mc = 2M^^^^ ■ Fc, (3) 

the factor 2 in the last equation represents two identical contributions Ci,C2 of the topology 
C. The amplitude then reads 

M = Mgorn " (^A + Fb + 2Fc) . (4) 

We turn now to the calculation of the form factors Fa, Fb, Fe- 
lt is convenient to use the standard Sudakov technique for calculating DL contributions, 
as described in [0, |^, and introduce Sudakov parameterization [|1^, [1^. One starts with 
decomposing the soft momenta in terms of those along the hard external momenta, for example 
ki, ^2 and transverse to them 

k = aki + [3k2 + k^. (5) 
The DL contribution comes only from the region 

mM/c_Lp < < s. (6) 

For different topologies it is convenient to use different decompositions of the soft momenta: 

2 9 

Tfl Tfl 

topology A: k = (3{j) p) + a{j) p) + k^ (7) 

s s 

topology B: k = [3k2 + aki + k± (8) 

2 

TfL 

topology C: k = [3k2 + a{p /C2) + k±_ (9) 

The loop integration in terms of the new variables reads 

00 00 00 00 

J d^k = ^ J da J d(3 J ndll. (10) 

—00 —00 —00 

The integration over the transverse momenta of the soft quark is performed by taking half of 
the residues in the corresponding propagator, 

f d'^k f s dadpd'^k_ 



1 1 /p xt^-n -^ -^"1 / d^d(3e{sap - m (11) 
J 2 sap — k i — X'' + lO 2 J 



k^-X^ + iO J 2sa(3-kl- X^ + iO 
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here A is a small (fictitious) gluon mass. In this manner the one loop amplitude may be 
calculated in the DL approximation and topology A gives, 

= / / ^e(a^ - ^)e{a - ^P)eiP - ^a) (12) 



CpOis / I , 2 1 1 



(-In^ — + ln — In— ), 13 



27r ^2 s s m? 

Topology B and C have following form at one loop in the DL approximation 



pi loop ^ pi loop ^ _C^jjdadp _rn^ C^A^^.r^^ 

27r J J al3 s 27r ^2 s ^ 





These results were derived first in [fzll and will be used in our calculation later on. 



2.2 The double logarithmic resummation 

The main idea behind our method is that the only origin of the double logarithms is the off- 
shell Sudakov form factor [] included in the effective one loop triangle diagrams A,B and C. 
Contributions from other types of diagrams will either cancel in subgroups or develop a simple 
on-shell Sudakov exponent due to the final quark anti-quark pair. The cancellations amongst 
the diagrams will be discussed in the next section where the general mechanism is elucidated. 
Here we study the DL resummation in the topologies A, B and C. 

The diagrams of topology A represent a simple case. Here we do not have any "hard" 
DLs, but only the "soft" ones, which develop the usual on-shell Sudakov form factor from the 
gluon exchanges between the final quark antiquark pair. It is well known, that the on-shell DL 
Sudakov form factor in QCD is the exponent of the one loop result 

Fa = Exp(f\-^"^), (15) 



with JP^^°°P given in eq.(|T2|). 

The resummation of the DL terms coming from the diagrams of topology B is more involved. 
First, we have to account for the "soft" on-shell exponent from the final quark antiquark re- 
scattering, similar to the one in topology A, eq. ([T5|) . Second, we observe a new element, namely 
the "hard" double logarithms from quark anti-quark re-scattering inside the one loop diagram. 
The interacting particles (quarks in this case) are hard and slightly off-shell, 2 ^ ■5- 

According to our general idea, the resummation of the "hard" DL can be obtained by 
including the off-shell Sudakov form factor into the loop, and accounting only for those momenta 
which do not destroy the DL nature of the one loop result. 

We recall that the off-shell Sudakov form factor, is a vertex of the production of the quark 



and anti-quark with small virtualities pf and P2 and has a form |13 



S{p„p,) = Exp( - ^ In(^) ln(^)), (16) 

assuming that ^ \P2\'^ s. Now the momenta pi and p2 become loop momenta, 

i.e., they depend on a,P 

pI = ih + = sf3, pI = {k2 -k) = -sa. (17) 
^ hard outer lines are slightly off-shell 



5 



Because of the DL approximation, the kinematic region of interest is restricted by the inequal- 
ities 

< bilMp2p < s. (18) 

Taking into account the result eq.(|T^) and including it into eq.(0) we get 

Fs = j j'-^e(a.a-'^)E.p{-^ln\a\inm) (19) 

Ztt J J ap s ^ Zn ' 



We next transform the exponent into a power series and find that the integral of the n— th term 
will be of the form 

/ ii. I rf&ff«er' = r(n + . + i)r(n + i + i) 

J ^1 y ^2^1 ^2 Y 3 + 2n + a + b) ^ ' 



i-a 

The final result at DL accuracy reads 



^ ^or(2n + 3)V ^ ^ 

with pb = ^2^" ' ~ (^)- "^^^ index n shows that the order of the amplitude is a". We 
can clearly identify the separate contributions of the fixed orders in as- On the other hand, 
if Pb is large all terms in the series are important, giving altogether some analytic function 
Fdl{p)- This function is identified with a hyper-geometric function 2-^2(1, 1; 2, |; 2;), namely 

Fb = i^r°°^E^|^(-/^B)"=2F2(l,l;2,|-^)Fr°- (22) 

These functions are exact answers in pb in the DL approximation. For large values of the 
parameter p the function F^l has the following asymptotic, 

Fb{pb) = (23) 
Pb 

Thus we see that despite the fact that perturbation theory blows up at large pb, the resummed 
result gives a smooth well defined function. 

Topology C differs from B only through the color structure at the DL level. At the end of 
the calculations the answers for topologies B and C are similar and are related by the simple 
substitution: 

Ci? — > Ca/2 in the variable pb- (24) 

Combining all results for the three topologies we have 

M = MBor^(l + Fb + 2Fc.)Exp(F^-'°°P). (25) 

with the functions defined above. We must stress that there are other diagrams which can give 
DL contributions. Fortunately, all of them cancel in the certain groups, as will be discussed in 
section 3. This cancellation has a very simple physical interpretation (see Sec. 3). 

Our results are written in two different representations. We have checked that they are in 
agreement with where the authors used a completely different method, not mentioning 
the off-shell Sudakov form factors at all. We believe that our approach is more transparent. It 
shows that the new double logarithms are of Sudakov type, and, therefore, the resummation 
procedure becomes simple. Another advantage of our approach is that it opens up the possibility 
for resumming the single logarithms as well. 
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Figure 2: A diagram responsible for determining the normalization scale of the coupling con- 
stant. 



2.3 Next-to-leading-logarithmic accuracy 

It is possible to develop this approach to achieve next-to-leading-logarithmic accuracy. 

First, the factorization formula should be modified in order to take into account single 
logarithms. The amplitude reads at NLL approximation 

M = Mbo^^(1 + A + Fb + 2Fc)Fa (26) 

with the function 

ln(s/m ), 



vr 2 

which can be extracted from the explicit calculations performed in by expanding the result 
at s = \t\ = \u\ 3> m^. Fa is the on-shell Sudakov form factor to the next-to-leading logarithmic 
accuracy 



Tl 



Second, in order to calculate the form factors Fb and Fc at NLL, we need an expression 
for the Sudakov form factor also to this accuracy. We turn now to the calculation of the form 



factors Fb and Fc- In fact, such an analysis already exists in the literature [|T5| for the case 
when the two external fermion lines are off-shell by the same amount, i.e., pf = P2 = P^'- 

SNLLip,p) = Exp( ^^^1^ (^) + ^""^J^^J- 

For our purposes, we need to extend the analysis to take into account that pf p\. It is easy 
to see that for the region, eq.(|l3), the proof of factorization and exponentiation given in 



goes through with straightforward changes. The major change involves the normalization of 
the coupling. We have studied one and two loop diagrams for the qq production vertex with 
slightly off-shell quarks with the following result: 

SNNLiPi,P2) = Exp[ — (ln^ln^ + -ln^ + -ln^Jj, (28) 



with the normalization of the coupling constant determined to be = ylPiUpil- We show 
only the double and single IR logarithms in the exponent of eq.f^Sf). 

In order to understand this normalization, we have to consider the diagram shown in Fig.0 
where we keep track of the nj dependent pieces only since they are separately gauge invari- 
ant. Such diagrams can be accounted for by considering the following gluon propagator 
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where n(/c^) is the vacuum polarization by the gluon; at the one loop level it is simply 11 = 
In (^^e*"^, /3o = 11 — \nf, C being a scheme-dependent constant {MS scheme C = — |). 

The diagram, Fig.^ corresponds to the first term in the expansion of the gluon propagator 
in as- The n^-part of this result is as mentioned earlier a gauge invariant subset of the complete 
set of two loop diagrams. 

Because the effects of the running coupling gives only single logarithmic terms it is enough 
to consider the remaining integrals to DL accuracy. Namely, we may trace only the terms 
proportional to In ^-^^ In from Fig.^ At the DL accuracy the spin structure of the amplitude 
is simple, so that one needs to consider the scalar integral only. 

To evaluate this, we consider a slightly more general integral 

j_ - f d'k 1 1 ^^^^ 

^ J (27r)4 (pi + ky {p2 - k) {k^y+^ ' ^ ' 

which after expanding in A will give us the desired integral I. Using Feynman parameters, this 
integral is reduced to 
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J = --^ I dy^fx''' e-AMA+i.) _ g-Ain(i.) ^ (32) 



Aln{A+B) _ -Aln(B) 



with 



{Airy J ^AA' 



A{y) = ply^ + 2p,p2y + pI (33) 
B{y) = -2pip2y - pjy - pi, (34) 
A{y) + B{y)=ply{-l + y). (35) 



The function A{y) has two zeros, y±: 



.21 V ^ -2pip2±^{2p^P2Y -ApIpI 

2 

For very small virtualities, p\-,p\ — > 0, the roots are simplified to y+ = — ^,?/_ = Ex- 
panding the integrand of J in A up to second order we have 

J = T^/fdy—^ \? r/x^^(in( ,^ plyi^-y) \ (37) 

(4vr)2i pliy - y+){y ~ y-r ^ ^{2piP2 + pl)y + pI^ 

+ ^[-\n'{ply{l-y)) + \n''ii2p,p2+pl)y + pl)]). (38) 
The final integration over y is simple, the result is 



J = y^-^^^^— -( -ln^ln^ + -ln^ln^ln(^^i^)). (39) 



(47r)22pip2 ^ s s 2 s s jj: 
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We see, that the first term in this equation reproduces the DL result from eq.(|l6D and eq . (^8]) . 
It can be checked, that the second term suggests the normalization of the couphng constant to 



be, = ylPillpil- Indeed, returning to the integral /, we find 



2tt 



In In 



It is clear that the last logarithm, containing the (3q term, can be absorbed into the running 



coupling, giving ^^(z/^) with the normalization point = y IpfHpil- The exponentiation of the 
integral / will give us the final off-shell Sudakov exponent, eq. (|28|) . 

In order to get single logarithms in eq. (p8D we have to include the numerator and the spin 
structure. We do not present the details of these calculations here. Instead we note, that all 
logarithms we have accounted for are of infrared origin, s Pi,P2 0. We do not show the 
UV logarithms which come as a result of the renormalization of the quark mass. Such terms can 
be omitted if the quark mass in the leading order result is normalized at a large scale //^ = s. 

The formula eq.(^8D reproduces the expression for the Sudakov form factor at non-equal 
virtualities at DL accuracy derived by Carrazone et. al. in ||13[, eq. (p!6D , as well as at NLL with 
equal virtualities p"^ = pi = pi derived by Smilga in JlSl , eq. 



In addition, the normalization point z/^ = y Ipil that we find reproduces that of the NLL 
results with equal virtualities p"^ = p\ = pi derived by Smilga in ||15[, eq. (pT]) . 

This scale, u"^ = \J\pi\\p2\, has a very transparent origin. The vertex of the interaction 
of a soft gluon with an off-shell quark (pf) is described by the coupling g{pl). In the situ- 
ation of gluon-exchange between two quarks with different virtualities, we have an effective 
coupling g{pDg{pl)- Using the running of the coupling g'^{fi'^) = 4:7ias{fi^), at one loop level, 
Q;(/i^) = a<i(z/^)/(l + ln(^)), we will find that the effective coupling g{pf)g{pl) is reduced 



to |pf 1 1^21), which coincides with our previous results. 

As a next step, we include this form factor inside the one loop triangle diagram and calculate 
the last loop integration with the form factor which now accounts for all large logarithms to 
NLL accuracy. The final result for the next-to-leading-logarithmic form factor reads 

F« = FE, + (41) 

with F^L from eq.(^T]) and 

_ 1 pi-ioop ^ r(n + l) ^ .nfo PbPo n ^ n + 1 ln(g//i^) ^^ 
^^^^ - L^^ „^or(2n + 2)^~^^^ 2n + 2 V2(2n + 3) + ^^^^ 



with Po = 11 — rif is the number of light flavors . 



Topology C at NLL gives a slightly different result. In the previous section where we have 
already studied the DL result, we saw that the "hard" Sudakov form factor is developed as a 
result of the re-scattering of a hard quark on the hard gluon. Let us start with the Sudakov 
form factor for the quark-gluon vertex with an off-shell outgoing quark of virtuality pf and a 
gluon of virtuality p^. The result for the reduced form factor is 

= £,p(-i!i(^(,„MtoM + bM + ii,M) (43) 

^ Atc ^ s s s 2 s ^ 

+ In — m 

47r ^ s ACf s 
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where the coefficient ai = — ^'^^"^ is related to the wave function renormahzation of the 
gluon. Now we can obtain a resummed result by substituting this form factor in the known one 
loop integral of topology C. The DL result coincides with that of topology B (as we discussed 
previously, after accounting for the substitution Cp ^ Ca/2). The result for next-to- leading- 
logarithmic form factor reads 

pC _ pi-ioop^ r(n+ 1) / 2Cf ai 

^^^^ - r(2n + 2) y^--Cl^2C; ^^^^ 

2pc/5o n f n + l ^ ln(s//i2) 



Ca 2ri + 2V2(2ri + 3) L 

The final result for the amplitude is eq.(^) together with eq.(^,|^,0). That concludes our 
derivation of the NLL form factors. 

We emphasize that throughout this subsection we have taken into account all logarithms 
of the Sudakov type only. In section 4 we argue that to the NLL level this is justified because 
of cancellations between diagrams that is similar to the DL case. In other words, logarithms 
of the NLL type coming from other sources (than the Sudakov form factors) cancel between 
contributions of subsets of diagrams. 



3 Cancellation mechanisms at the leading logarithmic 
level 

The simple expressions for the resummed results given in the previous section are a consequence 
of the cancellation of double logarithmic terms which are not of the Sudakov type. In fact for 
the exponentiation to work there has to be a large number of cancellations between subgroups 
of diagrams. Such cancellations have been discussed in the literature explicitly at the two loop 
[0 and the three loop levels |T^. We agree with the results of these papers however the new 



element we would like to add here is to point out the general mechanism which is responsible for 
these cancellations. We discuss this mechanism on various examples and write down identities 
which guarantee the cancellations to higher loop orders. 



3.1 The dipole mechanism for cancellation 

The relevant cancellations at the double logarithmic level are due to a single mechanism - the 
dipole mechanism. Two fermions with opposite charges form a dipole when separated by a 
short distance. When observing from far away, or probing by means of a gluon with long 
wavelength, one finds the dipole system to be neutral to leading order. To be more specific, 
the total amplitude for exchanging a soft gluon with two such oppositely charged fermions, as 
shown by the sum of the two generic Feynman diagrams in Fig. (^, vanishes in the infrared 
sensitive region (double and single logarithms). This infrared sensitive region is that of a soft 
gluon, i.e., a gluon with all components of the 4-momentum small, /c^ = (/c"*", /c_l), 

k+,k",k^_ ~ A ^ 0, 

This soft region is responsible for the double logarithms. 

The cancellation can be shown explicitly by employing the Grammar- Yennie decompo- 
sition and the Ward identity. To illustrate the mechanism let us consider a special case shown 
in. Fig. (^). The amplitude for the ffist diagram is 
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Figure 3: The general mechanism of the DL cancellation in the groups 1 and 2 of section 3.2 is 
the dipole interaction of a coUinear pair of quarks with one soft gluon. 



Mi = ri 



1 



-7 



1 



r. 



1 



{fi- ]i) - m (^1 - / - ^) - m (^1 + ^) - m 



(45) 



where we have explicitly written down only the factors relevant for the cancellation and denoted 
the rest of the corresponding expressions for the diagram by Fi and Their explicit form are 
of no interest as long as they are common to both the diagrams, Mi and M2. Employing the 
Grammar- Yennie decomposition we write for the gauge boson propagator, 



Daf3{k) = Gaf3 + K^f^ik), 

k^Dp^pl 



Gaf3{k) — Dap 



k-pi 

kgD^rPl 

k-pi 



(46) 
(47) 

(48) 
(49) 

(50) 



and keep only the term Kai3{k), i.e., we make the replacement 

n - So'f^ kag/3rPl 

which leaves the infrared behavior unchanged |]T^, we obtain 

{ki — ly — {ki ~ I — ky — m? (^1 + ^) — m /c^(pi ■ k) 

huh 



kaQprPl 



{ki — ly — mn? {ki — I — ky — nrfl 

2^1 (pi ■ ki) 

pi - ly - m2][(A;i -l-ky- m?] 



^2u{pi] 
^2u{pi) 



k'^ipi ■ k) 
1 

k'^ipi ■ k) ' 



(52) 
(53) 



where we have applied the Ward identity on the outgoing fermion line. Summing the two 
diagrams together, we have 



Ml + M2 = M 



Pi -ki p2- k2 



Pi-k p2- k 



(54) 



where the relative minus sign arises from the fact that the soft gluon interacts with two oppo- 
sitely charged fermions and 



M 



-Fi 



h 
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=0 



k2 P2 k2 

Figure 4: Example of a cancellation dictated by the dipole mechanism. 



In order to prove the sum, Mi + M2, vanishes, we use Sudakov parameterization, 

k = (3ki + ak2 + 
and proceed with the soft approximation, 

sa, s(3 ^ k^. 

We arrive at 



Pi ■ kl _ p2 ■ k2 
Pi ■ k P2- k 



1 



{f3 + a)^ - pi± ■ k± {(3 + a)^ - p2± ■ k_i 



(55) 



(56) 



(57) 



Again, in the soft approximation, pi± ■ k± and p2± ■ k± are negligible compared to (/3 + a) | due 
to the 5-function 

6{sal3 -kl- X^) (58) 
arising from the pole of the gauge boson propagator. This can be readily seen in the expansion 



1 



{(3 + a)^ - pi± ■ k± (/? + «) 



1 + 



Pi± ■ fcj 



+ 



Pi± ■ k± 

(/5 + «)f 



+ 



(59) 



since the odd terms in pi± ■ k± vanish when integrating over the angle between pi± and fcj 
Further, 



' 2"(f)"^Jcos2"^ a>/3 



Pi± ■ fcj 

(/? + «) 



2n 



2" 



a/3 



(a+/3)2 

n .,,2n 



cos^" e (3 



(60) 



. 2" (|)"^^cos2"^ /? > a 



and power suppression arises from the term {j^ with a ^ (3, with /? a and 

for all the three cases. 

Thus we have seen that the two diagrams indeed cancel each other in the double logarithmic 
approximation. The explicit calculation is applicable to the generic diagrams in Fig. (^ if we 
recall that what really matters here is nothing but the exchange of the soft gluon with the 
fermion-antifermion pair. The same conclusion holds for a exchanging coUinear gluon in the 
single logarithm approximation as we will discuss later. 
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3.2 Three loop examples 

We will next discuss how the dipole mechanism for cancellation works in certain 3 loop examples. 
More specifically, we will show the cancellation of double logarithms in the diagrams Zi — Ziq 
of topology C. For simplicity of presentation we first discuss only the Abelian case omitting 
all group theory factors and in the next section the extension to the non- Abelian case will be 
presented. These diagrams can be grouped in twos or threes according to the cancellation as 
seen below. 

1. Zi+Z2 = 

2. Z2 + Zi = 

3. Z5 + Zq + zj = 

4. zs + Z9 + zio = 

5. zn + Z12 + zi5 = 

6. Zi3 + Zu + ^16 = 

These six groups cover all the Abelian-like diagrams of Fig. |^. 
Group 1. 

We start with group 1. The amplitudes zi and Z3 may be written as 



Z3 



dkf 



dli 



(27r)4y (27r)4L ' ^^_/^_m^ + io' 



(62) 



where M^^ is some one loop sub-diagram, which is identical in both zi and 23. Such a rep- 
resentation of the diagrams is very convenient because it turns out that the cancellation of 
double logarithmic contributions could be traced without explicit computations of the separate 
diagrams. Instead we observe the cancellations by comparing diagrams in some subgroups and 
using kinematic simplifications specific to the double logarithmic approximation. 

Because the final quark line is hard we might write M^^^ = p2^M^. The momenta li must be 
softer than kj, otherwise no double logarithms can appear. Then we find that integrand of the 
diagram zi contains {p2 -kj). Therefore decomposing the momenta {kf) into kj and p2 we find 
that only the component parallel to kf contributes. We remind the reader that p\ = m?^ 0. 
Introducing n^, the unit vector parallel to kf, 



(63) 



we obtain for the sum of the diagrams 



Zl + Z3 



dkf 



dli 



(27r)4y (27r) 



■^M\n-p2) 



+ 



1 



nil + ^0 —nil + ^0 



= 



(64) 
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Figure 5: The Abelian diagrams of the topology C. 
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Thus the leading contribution in the integrand, responsible for the double logarithmic asymp- 
totic, cancels. This is shown in Fig. |^. 

As we will see below the mechanism of the cancellation is very general and is related to 
the dipole mechanism. Indeed we may define an impact factor as a sub-diagram, where a real 
photon produces a nearly on-shell quark and antiquark, which interacts with the soft gluons. 
The remarkable fact is that in the configuration responsible for double logarithms both particles 
move in the same direction. Because the quark momenta is harder than the gluon momenta, 
only the direction of its movement plays a role. 

kjp np 
kjl nl 

Further, because the hb system does not have a net color charge, being produced by a photon, 
the first term in the multi-pole expansion of interaction of this system with any number of soft 
gluons is the color-dipole moment. This leads to the power suppression of the sub-diagram and 
the loss of the double logarithms. 
Group 2. 

We continue with group 2 which constitutes diagrams Z2 and 2:4. The corresponding amplitudes 
could be written as 



Z2 



f ^^f 


f dli - 


1 {2nyJ 





(65) 



dkf 



dli 



(66) 



where as before, M^^^ is some sub-diagram, and we may write M^^ = p2^M'^. Again, the 
momenta li is softer than kj and only the component oikj parallel to kf contributes. We have 



Z2 + 2:4 



dki 



dli 



{2txY J (27r) 



P2) 



nil + ^0 —nil + iO 



= 



(67) 



We thus see that the leading contribution to the integrand, which gives rise to the double 
logarithmic asymptotic, cancels. 
Group 3. 

We turn now to the amplitudes for the diagrams z^, Zq, Zj which may be written as 



Z5 



dk 



f 



dli 



dU 



(27r)4y (27r)4y (27r)4 ^ + + ^ - + ^f+h-rriq + iQ 



-Z7 



Mfl (6J 



2^6 



dkf f dli f dl 



^7 



(27r)4y (27r)4y (27r)4 L^^^ + /i + - + iO ' ^j+Zi-m^ + iO 



Mfl (69) 



f ^^f 


f dli 


f dk 


1 (27r)4j 


(27r)4 J 


(2vr)4 



tY 



ll — TUq + iO ff + 1.2 — niq + ^0 



27 



(70) 
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Figure 6: The general mechanism of the DL cancellation in the groups 3-6 is the dipole inter- 
action of a coUinear pair of quarks with two soft gluons. 



it is clear from the diagrams that we can write to the accuracy needed for the relevant sub- 
graphs, Mj^l — Mj^l — Mj^l — p^ixPiuM^- The momenta li, are softer than kj and only the 
component oikj parallel to A;/ contributes. After simplifications, we obtain 



dl 



11 11 



(71) 



n{h + h) (nh) n{li + I2) {nh) {nh) {nh) 







Again, we see that the leading contribution to the integrand, which could give rise to the double 
logarithmic asymptotic, cancels. 
Group 4. 

The amplitudes for -Zsi -29, -^lo could be written as 



dk 



f 



dli 



dU 



?7 



(27r)4y (27r)4y (27r)4 L ' -m^ + iO ' - Ji^ - Ji2 - mq + 



(72) 



^9 



dkf 


f dli 


r dh 


(27r)4 } 


(27r)4 } 


{2nr 



h ~ h 



rrin 



^^K^ (73) 



ZlO 



f ^^f 


f dli 


f dk r 


1 (27r)4j 


(27r)4 J 


(27r)^ y 



M(;°) (74) 



with M(8) 



P2nPiuM'^ to the desired accuracy for the sub-diagrams. The 
momenta h, I2 are softer than kj and only the component of A;^ parallel iokf contributes. After 
simplifications, we obtain 



^8 



dU 



f dkf r dli f ^--2. n,.A, X/ X 

+ + = y (2^7 (2^0^/ (2^^'^^"-^^)^"-^^^ 



(75) 



1 1 



+ 



1 



1 



+ 



1 



{nh) {nh) n{h + h) {nh) n{h + h) {nh) 







We see that the leading contribution to the integrand cancels. The mechanism of the cancella- 
tion is identical to the one of group 3. It is important to note specially for purposes of the next 
section that = M^. Hence one might find another 2 sets of 3 diagrams in groups 3 and 4, 
which cancel in groups of three. 
Group 5. 

The amplitudes for zu, Z12, could be written as 



dh 



^ ( dkf_ f Jh_ r 
J {2n)^J {2%)^ J (27r)4 



Z7 



(76) 
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dk{ 



Zl2 



dli 



dlo 



(27r)4y (27r)4y (27r)4 L^^j + /i + ^ - + iO ' + + 



M(;,2) (77) 



dkf 



dli 



dU 



{27r)^J (27r)4i (27r)4 L^|^j + /i + ^ - + iO ' ^i-rUq + iO 



(78) 



with M^^J^ = = = PinPiuM^ the relevant sub-diagrams. The momenta li,l2 have 

to be softer than kf and only component of kf parallel to kf contributes. After simplifications, 
we obtain 



f dkf f dli f dl2 . 



dli f dh 



(27r)4y (27r)4y (27r)^ 



(79) 



1 1 



+ 



+ 



{nil) {nh) n{li + k) {nh) n{li + I2) {nl 







We again observe the cancellation of the leading contribution to the integrand. The mechanism 
of cancellation is identical to the one of group 3 and is the just the dipole mechanism. 
Group 6. 

Finally the amplitudes for ^13, ^14, ziq could be written as 



Z\Q — 



dk 



f 



dli 



{2t:Y J {2t:Y J (27r) 



dlo 



/,i-mq + iO }^f + /,2-'mq + iO 



Mlf (80) 



Zi4 



dkf 



dli 



dlo 



?7 



(27r)4y (27r)4y (27r)4 L ' + ' ^/ - /i - - + iO 



(81) 



Z12. 



dkf 



dli 



dlo 



?7' 



(27r)4 J (27r)4 J (27r)4 L ' ^^-mq + iO ' ^/ - /i - - + iO 



M(;,^) (82) 



where as before we write to this accuracy, Mj^^^ — M^^^ = M^^^ = pifj^pi^M^. The 
momenta li, I2 have to be softer than kf and only the component of kj parallel to kf contributes. 
After simplifications, we obtain 



Zl3 + Zi4 + ZiQ 



dk 



f 



dli 



dlo 



{2nrJ (2.)4^^'("-^^)^"-^^) 
11 11 

+ -T. TT^^ + 



(83) 



1 



1 



(nil) (n/2) n{li + I2) (n/2) n{li + I2) (nli) 







which explicitly shows the cancellation of the leading terms. The mechanism of the cancellation 
is identical to the one for the earlier groups and is the dipole mechanism. Note that for these 
two groups as well, = after interchanging li and l2- Hence one might expect to find 
another two sets of three diagrams in groups 5 and 6, which cancel in groups of three. 

In conclusion, we have shown that all diagrams zi — zie cancel in group of two or three due 
to the dipole mechanism. 
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3.3 Inclusion of non-abelian diagrams 

The inclusion of the non-abehan contributions do not pose any new difficulties. In fact, there are 
many simplifications because the final state quark and antiquark are produced by photons and 
hence carry no net color charge. When the color factors are included the cancellations take place 
between diagrams with the same such factors. For example, consider the cancellation between 
the diagrams of group 1. Here zi and have the same group factor since the only difference 
between them is an abelian vertex. The same is true for group 2. Thus the cancellation between 
the diagrams of groups 1 and 2 proceed as in the previous section even in the full non-abelian 
theory. 

Consider next, group 3 of the previous subsection. For this case we will discuss the cancella- 
tion in two different ways. In the first method we note that only the group theory factor for the 
diagram zq is different from that of Zr, and zy. Exphcitly, for zq the factor is TaTcTaT^TcTj, while 
for the other two it is TaTcThTaTcTf,. All other factors are the same as in the previous section 
with the rule that the group factors are overall multiplicative. In particular, the sub-diagrams 
Mj^l = = = P2fiPiuM'^ contain no color matrices. Now we can write the color factor 
of Zq as the one for z^ and Zj plus a left over term which is ifabcTaTcThiCp — Ca/'^)- Consider 
now the diagrams in group 4 of the previous section. Here the group factor associated with zs 
is TaT^UTaT^Tb which is different from zg and Zw which is T^T^Ta^T^T^. The difference now 
is —ifabcTaTcTbiCF — Ca/'2). Using the fact that = we see that the left over pieces 
from groups 3 and 4 cancel each other out. It is easy to check that the same applies to the 
combination of groups 5 and 6. 

It is clear from the above discussion, that it is much more straightforward in the non abelian 
case to change the assignment of the diagrams into different groups. It is easy to see that apart 
from the group theory factors, and in the soft approximation, the diagrammatic expressions 
for Zq and of the previous subsection (from groups 3 and 4) are the same (using of course 
= M*^). From groups 5 and 6 of the previous subsection the same applies to the expressions 
for Zi2 and Zi4^ (in this case using — M^). Thus the following assignment of the diagrams 
in the non-abelian case into different groups will make sure that all diagrams in a group have 
the same group theoretical factors: 

1. Zi + Z3 = 

2. z2 + z^ = Q 

3. Z5 + Zs + Z7 ^0 

4. Z6 + Z9 + Zio = 

5. Zu + Zu + 2^15 = 

6. 2^13 + Zi2 + ZiQ^O 

The new assignment docs not change the results for the abelian case and now the cancellation 
in the non-abelian case proceeds within each group. The only non-trivial result one must use 
is the equalities — and = M^. These are always seen to be true in the soft 
approximation. In the present grouping, the diagrams in each group are seen to be related to 
each other by a cyclic permutation of the gluon lines. Such a cyclic permutation leaves the 
color factor unchanged. 
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Figure 7: Non-Abelian diagrams. 
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Figure 8: The general mechanism of the DL cancellation is the dipole interaction of a collinear 
pair of quarks with many soft gluons. 



The cancellation between the other non Abelian diagrams shown in Fig.|^ also proceeds 
similarly: For example, Z21 cancels 2:22; -^23 cancels 2:24 and so on. 

This discussion now has been set up for generalization to all orders. 

3.4 Generalization to higher orders 

We see that the cancellation at the 3 loop level discussed in the previous subsection relies on 
the following: (1) the soft approximation, (2) algebraic identities like 

11 

+ z-ir = ^ (84) 



—p ■ ki p ■ ki 

11 11 11 

+ . . . . =0, 



-p ■ ki-p- {ki + k2) p- ki-p- k2 p- {ki + k2) p ■ k2 

and that diagrams in the same group are related by a cyclic permutation of the gluon lines. 
In the above, p is a generic hard momentum and ki are the soft ones. The soft approximation 
essentially tells us that a soft gluon does not see spin and more explicitly if p generically denotes 
a hard fermion momentum then we can consistently use, 2^ for the hard propagator and 2p^ 
for the vertex factors. This will ensure the equalities like = needed for the cancellations 
to occur for our process. As far as the algebraic identities are concerned, they are in fact special 



cases of the Sterman-Libby identities |T^. These identities obtained by considering diagrams 



related by cyclic permutations of the gluon lines read in general: 

n m—1 -I n -I 

E n — - — n — - — = (85) 

^0 t=0 P■(lm-P■q^ ,-4n+l P ' Qm ' P ' 

In the above, = /ci + ^2 + + ki and qi=o = 0. These identities are easily seen to reduce 

to Eqs. (|8^ for the cases n = 1,2. Thus we see each of the technical ingredients have a 
generalization to higher loop orders. The physics of the dipole mechanism and the color singlet 
nature of the final (and initial) states combined with the above guarantee the cancellations 
needed to all orders. 
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4 Contributing diagrams at the next-to-leading logarith- 
mic order 

The goal of this section is to identify those diagrams that either vanish or cancel some other di- 
agrams at next-to-leading logarithmic order. The diagrams left over are then just those needed 
in the discussion in Section 2.3. Our explicit analysis is only at the two loop level and at the end 
we argue that the results hold to NLL accuracy. In the first subsection, we discuss the regime 
contributing to the next-to-leading logarithms. We then introduce a power counting technique 
that enables us to discover the non- vanishing diagrams and make appropriate approximations. 
Armed with these two techniques, we are able to exclude a set of diagrams without actually 
carrying out the loop integrations. We remind the reader that we work throughout in the Feyn- 
man gauge. In this gauge we will see that for the process under consideration, the additional 
logarithms at the NLL level must be of coUinear origin. 



4.1 Sources of single logarithms 

In a two-particle scattering process, all momenta lie in the same plane. We can take two of the 
independent momenta, ki and /c2, as and '-' direction, and denote the third momentum as 
p {p could be either pi or P2) . The loop momentum / can be expressed in terms of light cone 
variables 

I^ah + Pk2 + 1±. (86) 
We note that the following integral in the soft regime 



/ 



dadp-^^e{saP - m^),i,j = -2, -1, 0, 1, 2, (87) 



gives terms hke log^ (^)"log ^ or (^)". In other words, the soft regime cannot give rise 
to a single logarithm at the one-loop level. 

Consider an n-loop Feynman integral with loop momenta k,! < i < n. We can decompose 
li in terms of 'external momenta' ki and /c2, 



h^aiki+ Pik2 + li±. (88) 



In the soft region, \ai\, ^ 1, the Feynman integral does not give the next-to-leading 

logarithm, log^""^ ^. This is because, intuitively, each integration over the light cone vari- 
ables a and (3 gives either a logarithm or a power suppression, as exemplified in the previous 
paragraph. We now prove this statement more rigorously. 

We first note the following replacement for loop momenta corresponding to a soft line, 

-5 \ : > -m5{sail3i + - m^). (89) 

tf — rrii + le 

The rest of the propagators can be categorized into hard, coUinear and soft ones. The hard 
propagators are irrelevant to the infrared sensitivity. The remaining possibilities are coUinear 
to ki or A;2, soft or coUinear to p. We discuss them in turn. 
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(1) Collinear to ki or k2. A boson that is parallel to ki has a propagator that can be written 



as 



1 



1 



2 1 



The fermionic propagator can be expressed as 

1 {h + Uili) - m 



{ki + UikY - ni^ 
h 



(2) Soft. A bosonic propagator 

1 



Using 6{saiPi + — mf), and 2 J (aiPj){ajPi) < ai(3j + aj[3i, it is easy to show 



(^Uiai)(^Uif3i)s\ > C^Uil 



Hence, the expansion below, 



{J2iUiai){J2i Uil3i)s 



1 



{J2iUiai){J2i Uil3i)s 



(90) 



(91) 



(92) 



(93) 



(94) 



converges. The fermionic propagator has a term proportional to /f , in addition to the expansion 
above. 

(3) Collinear to p. Again, the bosonic propagator can be expanded as 



{v + Y^iUikY - mf Ei^i("i + A) 1^1/2 + T,iUi{p± ■ h 



J:^u,ia, + f3i)\t\/2 



1 - 



E» Ui{p± ■ li±) 

EiM^(a^ + A)|t|/2 



(95) 



The fermionic propagator may give rise to an additional in the numerator. 

First, we only include the first term in the series in eq.(P^, and apply the following trick 



1 1 1 



1 



1 



aPa + P \a a + (3 I (3"^ 



(96) 



This trick reduces the number of different combinations of ct's and /3's while splitting one term 
into two. 

Therefore, when the propagators in case (1) and only the first terms in the expansions of 
the soft and coUinear-to-p propagators are included, the integral can be reduced to a finite sum 
of the type 



2n 1 
J soft regime • (J * 



(97) 
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where a's represent a's, /?'s or the combinations thereof. It is evident that an integration over 
each a gives either a logarithm or a power suppression m? / s. No next-to-leading logarithm can 
arise in the soft regime. 



We now include the whole series in eq.(0, |95|), as well as their fermionic counterpart when- 
ever appropriate. The denominators are polynomials of ai and j3i, 

111 

(98) 



Ei Uitti ' Y.i UiPi ' Ei Ui{ai + /3i) 
The numerators consists of terms proportional to p'^ and l^. The Feynman integral is thus of 
the following form 

f2{ai,/3i) 

where /i and /2 are both polynomial functions of their arguments and the spinor structure 
is not interesting. Since p± is the only vector in the integrand which is not the integration 
variable, the integral of interest can further be reduced to 

daidpid li± Y7 ^\ ' (100) 

J2(ai, Pi) 

where we have further left out the tensorial structure of pj_ and g^'^. 
By noting 

p±-k± = •\/piy^cos6'i (101) 

where the angles 6i are relative to the vector p±, the integral can be cast into 

I n dam E (102) 

where P and Q are polynomial functions of their arguments. We have used the 6 function 
(eq.(|5^) to perform the integration over and implicitly included the resulting 9 functions 
in the 'polynomials' P. The summation in the second line is due to the expansions in the soft 
and collinear-to-p propagators. It is a convergent series. 



Now we examine the arguments of the polynomial P in eq.( |102D . The \jY^{uiai + Vjj3j) 
represents various combinations of and (3j that may appear. For such a combination, we 
split the integration region 

/a-(l + i^ + ---), a, >/3, 



After such maniputions, we obtain a convergent series 

(104) 

with P and Q polynomial functions of ^/al, and In order to obtain the next-to-leading 

logarithm, 2?t, — 1 of the integrations have to give logarithm while the last one gives a constant 
of order 1. This is impossible for the above integral in the soft region. It follows that at least 
one of the loop momenta has to be taken out of soft region to correctly reproduce the next to 
leading logarithmic behavior in the Feynman gauge. 



23 



4.2 A power counting technique 



In order to identify the regime contributing to the next-to-leading logarithmic order at two- 
loop level, we consider inserting a gluon into the one-loop box diagram. From the previous 
subsection, we know the inserted gluon has to be collinear. Therefore, we will consider in turn 
the three topologies with an additional collinear gluon inserted to each of them. 

Throughout the following subsections, we always denote the soft loop momentum by / and 
that of the collinear gluon by k. 

We take two momenta p and p as the basis to decompose the momentum k of the collinear 
gluon, 

k = ap + (3p + k_i. (105) 

The generic momenta p and p can be any two of the external momenta ki, k2, pi and p2- 

In the so-called collinear region, without loss of generality, we assume k parallel to p such 
that 

lal ~ 1, 1/31 ~ W— < 1. (106) 
V s 

In general, all the propagators in a Feynman diagram can be characterized as hard (off- 
shell), soft or collinear (to a certain direction). In order to get the next-to- leading logarithm 
at two-loop level, which is a double logarithm multiplied with a single logarithm, a Feynman 
diagram has to contain at least four collinear and one soft propagators. Specifically, the double 
logarithmic form factor arises from a soft virtual particle 'interacting' with two (nearly) on-shell 
particles that are flying apart along two different directions. On the other hand, the interaction 
between two collinearly flying virtual particles gives rise to the single logarithmic form factor. 

As a result, whenever we have less than four collinear propagators, we can immediately 
conclude that the Feynman diagram does not contribute at the next-to- leading logarithmic level. 
When we have exactly four collinear and a soft propagator, we only keep terms proportional 
to a in the numerator. The (3 term can be dropped because it cancels a /3 in the denominator 
and effectively 'removes' a collinear propagator. An example of such a collinear propagator is 

1 1 . s 

(107) 



{p + k)"^ — m? (1 + a)f3{2p ■ p) + k'^ + p"^ + p"^ — m? 

And when we carry out the integration over fc, we can pick up a pole from an 'on-shell' prop- 
agator as the above one. It follows that f3s ~ k^ and k\ should also be dropped. If we get 
more than four collinear propagators, we will keep all the terms. However, in general the terms 
proportional to a are suppressed by rn? and thus leave j3 and k^ the leading terms. 



4.3 Vertex functions 

When inserting a gluon into a box diagram, we will obtain one-loop vertex subdiagrams inside 
four of the resulting two-loop diagrams. Two of the vertex corrections each have two legs 
(nearly) on-shell and the other soft. There is no large scale of order s, other than that from the 
UV cutoff, in such a sub-diagram. Therefore, these two vertex functions contain no infrared 
logarithms. 

The other two vertex function subdiagrams have two (nearly) on-shell and one off-shell legs 
each. They do contribute single logarithms in all the three topologies. Such vertex corrections 
are shown in Fig.^(g, h), |ll](g, h), [l^(e, f) and |13|(e, f). Self energy corrections are understood 
and not explicitly drawn. 
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ki pi ki pi ki pi 




Figure 9: Diagrams relevant to the next-to-leading logarithms in Topology A. The diagrams 
dependent only on Cp are shown here. 

Hereafter, we will omit the diagrams (and regions of diagram) that give rise to large loga- 
rithms only of ultraviolet origin, until we are ready to run the relevant parameters using the 
RGB. 

4.4 Contributing diagrams in topology A 

We show some of the diagrams relevant to the next-to-leading logarithms of Topology A in 
Fig.y. The shorter fermion lines in the box subdiagrams represent the off-shell propagators, 
which characterize the Topology A. Note Fig.|^(a) and (b) represent different regions of the 
same diagram, where the soft gluons are labeled by / and the coUinear gluons by k. The same 
comments apply to Fig]^(c) and (d). (We will follow the conventions that the characteristic 
off-shell propagator is denoted by the shorter line, / for the soft momentum and k for collinear 
momentum throughout the rest of this article.) 

The reduced diagram for the first six diagrams, Fig.|^(a) through (f), all consist of a hard 
vertex with four jets attached to it, as shown in Fig.|TD|. The two jets eventually emerge 
as quark and antiquark further interact with each other via a soft diagram (the gluon with 
momentum /). One of them consists of the collinear gluon and the quark or antiquark. In 
addition, these subset of diagrams are gauge invariant, since the reciprocal subset (consisting 
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Figure 10: Jet diagram for topology A 



of vertex correction and self-energy subdiagrams) are gauge invariant. We can now invoke 
results from a general power counting analysis of infrared sensitive contributions (both soft 



and collinear) to a typical wide-angle scattering process |T^, It was shown there that the 



logarithmic configuration requires that jet lines are attached to hard vertices by a single line, 



otherwise there is power suppression. The analysis of ||I9|, was made in a physical gauge, 
but it obviously holds for the gauge invariant set discussed above. Hence, the sum of Fig.^(a) 
through (f) does not contain infrared logarithms. 

The remaining diagrams are Fig.^(b) and (h), which contribute to the next-to- leading log- 
arithmic order. These are just the diagrams proportional only to Cp contributing to on-shell 



Sudakov form factor to NLL accuracy which is extensively discussed in fT^ 



4.5 Contributing diagrams in topology B 

Now we turn to the diagrams in Fig. ([lID . In Fig.([TT|e), the collinear gluon, labeled by k, can 
be parallel to pi, ki or k2. We discuss them in turn. 

i) Parallel to pi {k\\pi). There are only three collinear propagators left. This region can be 
excluded by power counting. 

ii) Parallel to ki {k\\ki). The soft fermion is the one labeled by /, while the fermion labeled 
hj k + I is collinear to ki. And the fermion with momentum k + pi is off-shell. There are 
exactly four collinear and one soft propagator left. Hence we only keep the component of k 
that is parallel to ki. The numerator of the diagram here is, 

N oc {h + ^ + Mh)i^ + ^) 

oc hKh)h = 0, (108) 

which implies a vanishing contribution to the next-to-leading logarithmic order from this region. 

iii) The last possible region is fc||/c2, which vanishes due to the similar reason as in ii). 
Therefore, Fig.|ll](e) vanishes as well as Fig.|Tl|(f). Diagrams Fig.|Tl|(a) - (d) can be shown 

to factorize. Take Fig.|Tl](d) as an example, the numerator of the amplitude is 

Nied = Hpih''{h + ^ + li^)Y{h + MkiW{k2){-h + ^ha{-h + lthAP2) 

f n{Pih''ih + MkiMik2){-h + fha[-^P2-{ki + k)]v{p2), k\\k^ 

\ U{pi)r [4fcl ■ {-P2 + k)] ih + MhW{k2){-h + ^)laV{p2), k\\p2 ^ ' 

The factorization is evident now. Similar results hold for the other three diagrams. However, 
they do not contain any large single logarithms. 

The remaining diagrams, Fig.|TT|(g) and (h) are just the ones included in the resummation 
discussed in Section 2.3 to this order. 
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Figure 11: Diagrams relevant to the next-to-leading logarithms in Topology B. The t-channel 
diagrams dependent only on Cp are shown here. 
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Figure 12: Diagrams relevant to the next-to-leading logarithms in Topology C. The s-channel 
diagrams dependent only onCp are shown here. 



4.6 Contributing diagrams in topology C 

In topology C, the two-loop s-channel diagrams only proportional to Cp are shown in Fig.(|l2D. 
The regions contributing to the next-to- leading logarithmic approximation for Fig.(|T2|a, b, c) 
are both the gluons being parallel to p2- 

Note the 'incoming' quark and the 'outgoing' gluon of the hard subprocess 7g qg are 
nearly on shell. In addition, the gluon labeled by k in the three diagrams are nearly on-shell 
too. We can expect the contributions only proportional to Cf of the three diagrams to cancel 
to the leading order, in the same manner as discussed earlier in Section 3. 

In order to show the cancellation, we decompose 

k = akk2 + PkP2 + k±. (110) 
Each of the three amplitude takes on the form 

Mi = M.^ ,^ ,, iAP2) (111) 



''2 



= -"^^ — {^;tw — 

with Mi, i = 1,2,3, representing amplitudes of Fig.(^a, b, c) respectively. In this subsection, 
integrations over / and k should be understood in the amplitudes. We consistently omit common 
numerical factors for simplicity. To the single logarithmic approximation, and k_i_ terms in 
the numerators can be neglected. Hence, 

Pk {p2 + ky 
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Figure 13: Diagrams relevant to the next-to-leading logarithms in Topology C. The u-channel 
diagrams dependent only on are shown here. 



and we have put back the afc/c2 + k± in the second line. 

The summation of k^Mi^ in the three amplitudes closely follows the earlier proof using the 
Ward Identity (see section 3.1). We obtain (for the piece proportional only to Cp) 



M1 + M2 + Ms 



■/(^2)7c;r7A^^(p2 



1 



1 



(113) 



The first term is suppressed by the quark mass, while the second one vanishes in the single 
logarithmic approximation by simple power counting. Therefore, the contributing diagram in 
Fig.^ are d, e, and f. 

Similarly, the diagrams with the u-channel hard subprocess also include six diagrams. Note- 
worthy is that Fig.p!3|(a) can also be drawn as a vertex correction to the s-channel diagrams, 
with the triangle subdiagram being the vertex correction therein. However, we notice they 
represent different regions. In Fig.|l3|(a), the quark labeled by / is soft, whereas in the s-channel 
diagram, the corresponding one is collinear to /c2- A similar remark applies to Fig.|T3](b). Here, 
it can be viewed as a vertex correction to the u-channel diagram, which in turn represents a 
different region. As a summary, the two diagrams shown in Fig.|T^(a, b) represent the region 
that both the virtual gluons are parallel to p2, while in their counterparts, the two virtual 
gluons are parallel to ^2 and p2 respectively. Obviously, the sum of Fig.O(a, b, c) vanishes to 
the next-to-leading logarithmic approximation in exactly the same way Fig.0(a, b, c) do. 

Therefore, the contributing diagrams of the Topology C at the two-loop level include 
Fig.0(d, e, f) and Fig.|T3|(d, e, f) only. We note the diagrams contributing to the "hard" 
off-shell Sudakov form factor in this topology are mostly non-abelian which are of Sudakov 
type. Some of these are shown in Fig.|T^. This completes our argument justifying that the only 
source of logarithm at this order are of the "Sudakov" type. 
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Figure 14: Some of the diagrams proportional to Ca in topology C. 



4.7 Extension to higher loops 

In this section, we argue without detailed proof that the above results hold to the full NLL 
accuracy. We consider only topologies B and C, since topology A, which corresponds to the 
on-shell Sudakov form factor, has been discussed in detail in |jl4 . 



For topologies B and C, consider the insertions of a gluon into the bare diagram. There are 
two cases: 

i) The gluon momentum is in the soft region. The analysis in Section 3 can be applied and it 
can be factored out. 

ii) If the gluon momentum is in collinear region, the analysis in the previous subsections apply 
and again, factorization results. 

All subsequent insertions of gluons in case (ii) must be restricted to the soft region. There- 
fore, for these the analysis of section 3 apply. In case (i), we keep inserting gluons and apply 
the same analysis until we encounter a collinear gluon. Then case (ii) applies. 

Hence we can conclude that to NLL accuracy, the relevant diagrams are all of the Sudakov 
type. 



5 Summary and Discussions 

In this paper we have studied the resummation up to the next to leading logarithmic level of 
the QCD radiative corrections to bb production by photon photon collisions. Apart from the 
phenomenological applications, this problem has inherent interest in providing a theoretical 
laboratory to study QCD effects. We showed that to the accuracy considered all logarithms 
are of the Sudakov type. On a diagram by diagram basis other types of diagrams do give rise 
to next to leading logarithms but they cancel amongst each other by the dipole mechanism. 
We explicitly showed how the dipole mechanism works to 3 loops and outlined an all orders 
generalization. At the NLL level we found that in the Feynman gauge the only new logarithms 
are of collinear origin. The final expression for the resummed amplitude is given in section 2, 
eq.(ll) together with eq.(El1,E2|,p[). 
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